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Observations with the Imaging Atmospheric Cherenkov Telescopes H.E.S.S., MAGIC, 
CANGAROO III and VERITAS have shown that the Universe is more transparent than 
expected to gamma rays above 100 GeV. As a natural explanation, the DARMA scenario 
has previously been proposed, wherein photons can oscillate into a new very light axion- 
Uke particle and vice- versa in the presence of cosmic magnetic fields. Here we demonstrate 
that the most recent observations further support the DARMA scenario, thereby making 
the existence of a very light axion-like particle more likely. 

1 Introduction 

A generic prediction of many extensions of the Standard Model of particle physics - including 
four-dimensional theories, compactified Kaluza-Klein models and superstring theories - is the 
existence of very light axion-like particles ( ALPs) . They are defined by the low-energy effective 
lagrangian 

Calp = Id^ad^a-^m'a'-^F^'^F^^a, (1) 
where F^" is the electromagnetic field strength, is its dual, a denotes the ALP field whereas 

]^/2 — 1/2 ~l/2 

TO stands for the ALP mass. Accordingly, it is assumed M ^ Gp and m <^ Gp {Gp 
250 GeV is the Fermi constant). The standard axion is the most well known example of ALP, 
but as far as generic ALPs are concerned the parameters M and to are to be regarded as 
independent. 

While extremely elusive in laboratory experiments, ALPs can give rise to dramatic astro- 
physical effects owing to the characteristic 77a vertex in >Calp ■ Correspondingly, ALPs can be 
emitted by astronomical objects of various kinds, and the present situation can be summarized 
as follows. The negative result of the CAST experiment designed to detect ALPs emitted by 
the Sun yields the bound M > 0.86 ■ 10^° GeV for to < 0.02 eV. Moreover, theoretical consid- 
erations concerning star cooling via ALP emission provide the generic bound M > 10^" GeV, 
which for to < 10~^^ eV gets replaced by the stronger one M > 10^^ GeV even if with a large 
uncertainty. The same 77a vertex produces an off-diagonal element in the mass matrix for the 
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photon- ALP system in the presence of an external magnetic field B. Therefore, the interaction 
cigcnstates differ from the propagation eigenstates and photon- ALP oscillations show up. 

A few years ago, it was realized that photon- ALP oscillations in cosmic magnetic fields - the 
so-called DARMA scenario [T] - can provide a natural explanation for the anomaluosly large 
transparency of the Universe above 100 GeV observed since 2006 by the Imaging Atmospheric 
Cherenkov Telescopes (lACTs) H.E.S.S., MAGIC, CANGAROO III, VERITAS, and further 
confirmed by subsequent observations. 

Our aim is to show that the most recent observations of blazars at rcdshift z > 0.3 beautifully 
fit within the DARMA scenario, thereby providing further support in favour of the existence 
of an ALP with M slightly larger than 10^^ GcV and m slightly smaller than 10~^°eV. Re- 
markably enough, this claim can be tested with some detectors originally devised to search for 
non-baryonic dark matter [2] as well as with future photon regeneration experiments [3]. 

2 Extragalactic background light 

Photons from distant sources scatter off background photons permeating the Universe, thereby 
disappearing into electron-positron pairs and so giving rise to a cosmic opacity. The corre- 
sponding cross section a{jj — >■ e^e~) peaks where the very high energy (VHE) photon energy 
E and the background photon energy e are related by e ~ (500 GeV/i?) eV. As far as lACT 
observations are concerned, the cosmic opacity is dominated by the interaction with diffuse 
background photons with 0.005 eV < e < 5eV, usually called extragalactic background light 
(EBL). Owing to the absorption process in question, photon propagation is controlled by the 
optical depth r(£', z), with z denoting the source redshift. Within the standard Big Bang model 
we have E[z) = Eo(l + z) and e(z) = eo(l -I- z), with Eq and eo referring to the present [z = 0). 
Therefore, the observed photon spectrum (f>obs(£'o, z) of a source at z is related to the emitted 
one $cm(£'(z)) by 

$obs(So, z) = e---(^'"^) (So(l + z)) . (2) 

Note that t^{Eq,z) increases monotonically with z, since a greater source distance entails a 
larger probability for a beam photon to be absorbed. 

Given the fact that all blazars observed so far by lACTs lie in the energy band 0.2 TeV < 
Eq < 2 TeV, we restrict our attention to this energy range from now on, in which blazar data 
are fitted as <f>obs oc E~^°'^'. The results for all VHE blazars observed to date by lACTs are 
exhibited in Fig. 1, where the corresponding observed spectral indices Fobs are plotted vs. the 
source redshift z for all sources reported in Table 1 of ref. [4] plus three more recently detected 
ones, namely S5 0716+714 at 2 0.31, 3C66A at z = 0.444 and PKS 1424+240 at z < 0.66. 

It is generally assumed that $cm oc E^^™ for E > 100 GcV, and so Eq. ([2]) yields 

F^P(z)(xFe,n + r^(So,2) , (3) 

up to a negligible additional term with a logarithmic z-dependence and z-independent factors, 
where CP means that this result follows from conventional physics alone. 

A key ingredient in the evaluation of Tj{Eq,z) is the EBL spectral number density, which 
depends on the adopted model for the EBL. We choose the one of Franceschini, Rodighiero and 
Vaccari (FRV) [5] . A convenient analytic fit to the corresponding EBL spectral number density 
at z = is 

n^(eo, 0) ~ 10-3 a (^) cm-^ eV^ , (4) 
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with 0.5 < a < 3. Besides redshifting all energies in proportion of 1 + z, the cosmic expansion 
dilutes the EBL by a factor {1 + z)^. In addition, the EEL spectral energy distribution changes 
because of the intrinsic evolution of the galactic population over cosmic times. A quantitative 
analysis shows that the EBL photon number density acquires an extra factor (1 + z)"^-"^ as 
long as z < 1. Altogether, we get n^(e(z), z) ~ (1 + z) ' n^(eo, 0). We expect the EBL photon 
absorption to be negligible for nearby blazars (z < 0.03), and so we suppose that observations 
of these sources do yield $om(£'). We find Fem — 2.4 on average. We further assume that all 
VHE observed blazars have an emission spectrum with basically the same slopB- Then the 
resulting values of r^,^(z) for 0.5 < a < 3 he between the two dotted lines in Fig. 1. 

As it is clear from Fig. 1, the actually observed spectral index Fobs(-z) increases more slowly 
than r^^(z) for redshifts z > 0.2. Moreover, the observed values cannot be explained for 
z > 0.3 by the EBL model of FRW even for a as low as 0.5. Being t^{Eo, z) a monotonically 
increasing function of z, we interpret the conflict between Fobs (-2) ^^nd F^|^(z) shown in Fig. 1 
as calling for a departure from the conventional view. 

A suggested way out of this difficulty relies upon the modification of the standard emis- 
sion mechanism. One option invokes strong relativistic shocks [5]. Another rests upon photon 
absorption inside the blazar [7] . While successful at substantially hardening the emission spec- 
trum, these attempts fail to explain why only for the most distant blazars does such a drastic 
departure from the conventional view show up. 



3 DARMA scenario and its predictions 

The spirit of the DARMA scenario [l] is quite different. Implicit in previous considerations is the 
hypothesis that photons propagate in the standard way throughout cosmological distances. We 
suppose instead that photons- ALP oscillations occur in the presence of cosmic magnetic fields, 
whose existence at the nanogauss level has been suggested by AUGER observations [8]. Once 
ALPs are produced close enough to the source, they travel unimpeded throughout the Universe 
and can convert back to photons before reaching the Earth. Since ALPs do not undergo EBL 
absorption, the effective photon optical depth gets reduced, thereby leading to a reduction of 
the predicted observed spectral index F|j'|^^'^^(z) in this context. 

Owing to the notorious lack of information about the morphology of cosmic magnetic fields, 
one usually supposes that they have a domain-like structure [9]. That is, B ought to be 
constant over a domain of size Ldom equal to its coherence length, with B randomly changing 
its direction from one domain to another but keeping approximately the same strength. As 
explained elsewhere [5], it looks plausible to assume the coherence length in the range 1 — 10 Mpc. 
Correspondingly, the inferred strength lies in the range 0.3 — 1.0 nG [5]. 

Following the same computational procedure as in ref. [1], we evaluate the probability 
P-y^-y(Eo, z) that a photon remains a photon after propagation from the source to us when 
allowance is made for photon- ALP oscillations as well as for photon absorption by the EBL. As 
a consequence, Eq. ([2]) becomes 

$obs(£^0, Z) - P-f^iiEo, Z) $em (£^o(l + z)) (5) 

so that Eq. ^ gets replaced by 

r^^(^) oc Fe„, - In P-,^,{Eo, z) , (6) 

^Note that no assumption is being made about the intensity of the VHE emission, so that we are absolutely 
not supposing VHE blazars to be standard candles. 
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again up to a negligible additional term with a logarithmic z-dependence and z-independent 
factors. Assuming m < 10"^° cV, M ~ 4 • lO" GcV, B = 0.5 nG and Ldom = 7Mpc as in rcf. 
[1] and adopting the same EBL model as before as well as Fcm — 2.4, the resulting values of 
r^j^^^^(z) for 0.5 < a < 3 lie between the two sohd lines in Fig. 1. 

Clearly, by a small change in the free parameters there is no problem to slightly lower the 
latter region so as to achieve a perfect agreement with observations. 

More details and additional references can be found in ref. [1]. 




Figure 1: The observed values of the spectral index for all blazars detected so far in the VHE 
band are represented by big dots and corresponding error bars. Superimposed on them is the 
predicted behaviour of the observed spectral index within two different scenarios. In the first 
one (area between the two dotted lines) F^^ is computed in terms of conventional physics in the 
FRV model of the EBL. In the DARMA scenario (area between the two sohd hues) F^^'^'^^ is 
evaluated within the proposed photon- ALP oscillation mechanism as based on the same FRV 
model of the EBL. 
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